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The existence of current-time universe’s acceleration is usually modeled by means of two main 
strategies. The first makes use of a dark energy barotropic fluid entering by hand the energy- 
momentum tensor of Einstein’s theory. The second lies on extending the Hilbert-Einstein action 
giving rise to the class of extended theories of gravity. In this work, we propose a third approach, 
derived as an intrinsic geometrical effect of space-time, which provides repulsive regions under certain 
circumstances. We demonstrate that the effects of repulsive gravity naturally emerge in the held 
of a homogeneous and isotropic universe. To this end, we use an invariant definition of repulsive 
gravity based upon the behavior of the curvature eigenvalues. Moreover, we show that repulsive 
gravity counterbalances the standard gravitational attraction influencing both late and early times 
of the universe evolution. This phenomenon leads to the present speed up and to the fast expansion 
due to the inflationary epoch. In so doing, we are able to unify both dark energy and inflation in 
a single scheme, showing that the universe changes its dynamics when -A = —2 H, at the repulsion 
onset time where this condition is satisfied. Further, we argue that the spatial scalar curvature 
can be taken as vanishing because it does not affect at all the emergence of repulsive gravity. We 
check the goodness of our approach through two cosmological fits involving the most recent union 
2.1 supernova compilation. 


Introduction General relativity predicts a universe 
characterized by a decelerated expansion, if only dark 
matter and baryons are involved within the energy- 
momentum tensor. Observations, however, predict a 
late-time accelerated universe, pushed up by some sort 
of exotic constituent whose nature is today unknown [1]. 
To account for this experimental fact, cosmologists em¬ 
ploy in general two different strategies. The first one 
consists in extending the Hilbert-Einstein action, mod¬ 
ifying general relativity by means of additional degrees 
of freedom Q|. The second one takes into account an 
additional fluid, dubbed dark energy , which enters Ein¬ 
stein’s equations and represents a source for the cosmic 
speed up [|j]. Both those treatments are based on the 
hypothesis that space-time geometry is fueled by compo¬ 
nents entering the energy-momentum tensor T^ v . In such 
a way, the energy-momentum tensor becomes the source 
for gravitation and modifies the geometry itself. How¬ 
ever, if we rewrite Einstein’s equations in an equivalent 
way as — \g^T = it is possible to reverse the 

interpretation of T^ u . Indeed, one can address a dynam¬ 
ical problem by considering an energy-momentum tensor 
induced by geometry, and not vice-versa. The philosophy 
is completely different from the one used in higher-order 
theories of gravity, since no extensions are required, but 
only a physical mechanism inside the Einstein equations 
themselves, providing a geometrical source term. 

It is therefore possible to imagine that under certain 
circumstances the geometry could correspond to a repul¬ 
sive field that counterbalances the attraction of gravity 
and induces an accelerated universe. This third way of 


handling the universe dynamics becomes a natural conse¬ 
quence of general relativity itself, instead of a new addi¬ 
tional ingredient that must be added into the energy- 
momentum tensor. Indications of repulsive effects in 
gravity are not new, specially in the case of naked singu¬ 
larities Q. Several intuitive approaches have been pro¬ 
posed to understand repulsive gravity, but only recently 
an invariant definition in terms of geometric objects was 
formulated [§]. Indeed, it was shown that the behavior of 
the curvature tensor eigenvalues can be used to find the 
places where repulsive gravity appears. This was explic¬ 
itly shown in the case of naked singularities with black 
hole counterparts where it was found that there exists 
a repulsion region which is always situated very close to 
the source of gravity. 

Here, we extend this approach to the case of a homoge¬ 
neous and isotropic universe, modeled by the Friedmann- 
Robertson-Walker (FRW) metric. We consider the effects 
of repulsive gravity as a source for the cosmic speed up, 
and find the stage at which the universe starts accelerat¬ 
ing, corresponding to a phase in which repulsive gravity 
dominates over attraction. We therefore interpret the ac¬ 
celeration of the universe in terms of a natural outcome 
due to the non-linearity of large-scale cosmology which is 
reflected in the eigenvalues of the curvature tensor. We 
characterize gravitational repulsion by evaluating precise 
values for the observable acceleration parameter q and 
its variation j. We also find that the universe dynamics 
does not depend upon the scalar curvature; interestingly, 
this could be interpreted as an explanation of the exper¬ 
imental evidence for that the universe appears spatially 


2 


flat at all stages of its evolution, and enables us to unify 
the inflationary phase with the accelerated phase in a 
single description. Finally, we find an approximate effec¬ 
tive density induced by repulsive gravity, assuming that 
the conditions valid at the moment of the repulsion on¬ 
set continue holding at all times. We check the goodness 
of our approach through a numerical test based on su¬ 
pernova union 2.1 data, performed on a grid in which 
all cosmic parameters evolve. Particularly, we employ 
two hierarchical experimental tests. In the first one, we 
take all coefficients free to span in allowed prior domains. 
In the second one, we fix the acceleration parameter by 
means of the most recent Planck results. We definitively 
show compatible limits between our theoretical predic¬ 
tions and numerical outcomes. 

Repulsive effects in a homogeneous and 
isotropic universe Although there are several intu¬ 
itive ways to define repulsive gravity, only recently an in¬ 
variant definition in terms of the curvature tensor eigen¬ 
values was formulated which gives physically reasonable 
results in the case of all naked singularities investigated 
so far Q. The idea is that the curvature tensor eigen¬ 
values contain in an invariant manner all the informa¬ 
tion about the behavior of gravity. To this end, it is 
convenient to introduce an orthonormal frame d a which 
represents the simplest choice for an observer to perform 
local time, space and gravity measurements so that all 
the quantities related to this frame are coordinate in¬ 
dependent. The orthonormal tetrad is determined by 
the relationships ds 2 = g^dx^dx 1 ' = rj a bd a ® d b , with 
r] a b = diag(—1,1,1,1), and d a = e a dx^. Furthermore, 
the first and second Cartan equations enables one to com¬ 
pute explicitly the Riemann curvature components in this 
frame. 

For the analysis of the eigenvalues it is convenient to 
consider the bivector representation of the Riemann ten¬ 
sor, which is also useful for finding its irreducible rep¬ 
resentation with respect to the Lorentz group. We use 
the notations and conventions introduced in |7[, accord¬ 
ing to which each bivector index A,B = 1,..., 6 cor¬ 
responds to two tetrad indices, i.e., A —► ab , with the 
convention 1 — > 01, 2 —> 02, 3 —> 03, 4 —> 23, 5 —> 
31, 6 —> 12. Then, the curvature tensor can be expressed 
as a (6 x 6)—matrix which, in the case of the FRW or¬ 
thonormal frame d° = dt , d 1 = a(t)dr/y/l — kr 2 , d 2 = 
affjrdO, d 3 = a{t)r sin ddip, reduces to the diagonal ma¬ 
trix 


n = diag(A, A, A, r, r, r) , (1) 



a 


- g (3-P + p) , 



( 2 ) 


According to our invariant definition of repulsive grav¬ 
ity [fjj, the curvature eigenvalues reflect the behavior of 
the gravitational interaction and if gravity becomes re¬ 
pulsive in some regions, the eigenvalues must change ac¬ 
cordingly; for instance, if repulsive gravity becomes domi¬ 
nant at a particular point, one would expect at that point 
a change in the sign of at least one eigenvalue. Moreover, 
if the gravitational field does not diverge at infinity, the 
eigenvalue must have an extremal at some point before 
it changes its sign. This means that the extremal of the 
eigenvalue can be interpreted as the onset of repulsion. 

Let us now consider the eigenvalues of the FRW cur¬ 
vature tensor. First, we note that r is positive definite 
by virtue of the energy conditions and so it does not con¬ 
tain any information about repulsive gravity. It therefore 
represents a constraint eigenvalue which fixes the total 
energy budget pushing up the universe, but with no in¬ 
formation on its evolution. 

Instead, A can change its sign in the interval where the 
pressure P is negative and, so, the universe can acceler¬ 
ate. Thus, A represents the dynamical eigenvalue which 
dominates the spacetime dynamics, predicting the accel¬ 
eration. Moreover, A does not depend explicitly on the 
value of the constant k which is the first interesting con¬ 
sequence of our approach. Indeed, this means that the 
universe does not need scalar curvature to accelerate so 
that the curvature density flk = — a k H 'i can be taken ar¬ 
bitrarily small, without affecting the universe dynamics. 
This theoretical result is in agreement with observations 
0- On the other hand, this could be considered as a 
theoretical explanation of the spatial flatness of the uni¬ 
verse. 

To further analyze the repulsive effects in the dynamics 
of our tire universe, we introduce the acceleration param¬ 
eter 


m 


a _ H 

aH 2 ~ H 2 ’ 


(3) 


to obtain 


A = qH 2 (4) 

so that A < 0 when the universe is accelerating, whereas 
for A > 0 the universe decelerates. The condition 

H = -H 2 (5) 

determines the special point A = 0 at which the acceler¬ 
ation parameter vanishes and the universe undergoes a 
transition from a decelerating phase into an accelerating 
phase. The particular time t = t acc at which this phase 
transition occurs follows from the integration of Eq. 
i.e., 


Hacc — H(t acc ) 


1 


where we have used the Friedmann equations and 87 tG = 
1. It then follows that the only independent eigenvalues 
of the curvature tensor are A and r. 


tacc T Cq 


( 6 ) 








3 



0.0 0.2 0.4 0.6 0.8 1.0 


Z 

FIG. 1: Plot of X in the redshift domain z G [0; 1] with H and 
q from the KCDM paradigm with = 0.318 and ho = 0.68. 

where Cq is an integration constant. The behavior of the 
eigenvalue A is illustrated in Fig. Q] for q and H of the 
ACDM model. 

As mentioned above, the onset of repulsion is de¬ 
termined by an extremal of the eigenvalue, i.e. A = 
qH 2 + 2qHH = 0. Then, using Eq. 10, we obtain the 
onset condition 


q = 2H(l + q)q. 

(7) 

Using the identity [ 13 ] 


Q = -(j -q-2q 2 )H , 

(8) 


where the jerk parameter j = 'a / (aH 3 ) determines the 
change of the acceleration [13], the repulsion onset con¬ 
dition © reduces to 

j = ~Q , (9) 

implying that at that moment of time t = t rep the jerk 
must be positive because the onset happens when the 
universe is still decelerating (q > 0). By using the def¬ 
initions of the cosmographic parameters, we can rewrite 
@ as 

f = -2 H, (10) 

a relationship that can be integrated and yields 

H r ep — H ( trep ) = a tanh[a (t rep + /?)] , (11) 

where a and /3 are integration constants. This is the be¬ 
havior of the Hubble parameter at the onset of repulsion. 
From that point on, the value of the acceleration param¬ 
eter q diminishes until it vanishes at t = t aC c (cf. Fig. [lj. 
It seems reasonable to construct a cosmological model 
based upon the behavior of the Hubble parameter near 
the repulsion onset point CD. where the acceleration be¬ 
gins to play an important role in the universe dynamics. 


In fact, one expects acceleration to be important in late 
and early phases of the universe evolution. Hence, we 
can assume that Eq. CD is valid from the transition time 
up to z —> 0 and also during the inflationary phase. 

As the redshift decreases, corresponding to the limit of 
t —> to, with to our present time, we get 

a tanh[a (to + /3)] ~ pa , (12a) 

a 2 costT 2 [a (t 0 + /?)] « -3 H 0 pm - Pm , (12b) 

where we considered pA = 3A the cosmological constant 
density in units in which 87rG = 1, p m is the today matter 
density and, furthermore, we used the continuity equa¬ 
tion p+3H(P+p) = 0 for the case of a late-time universe, 
in which P ~ Pa and p ~ p m + pa ■ The above equations 
fix limits over a and /? at late times. 

On the other hand, at the inflationary phase we have 

atanh[a (U + /3)] ss V(<j >), (13a) 

a 2 cosh -2 [a (U + /?)] « — \/3 V (</>)</>, (13b) 

where we considered the relation pa: p$ = 1 <f> + V((f>) 
and P a P^ = V((/)), with cf> and V((j>) the field and 

potential determining the inflation, respectively. We also 
took into account the slow-roll approximation V(4>) (f> 

and ti , the time at which inflation occurred. Those rela¬ 
tions impose bounds on a and /3 at the times in which the 
universe accelerated. In other words, the mechanism of 
acceleration is induced for different epochs in the same 
way, but with different coupling constants, whose val¬ 
ues are associated to the typology of the acceleration in¬ 
volved, i.e., with dark energy or inflation. 

Numerical outcomes from experimental fittings 
We carry out two analyses based on supernova standard 
candles, employing the most recent union 2.1 compila¬ 
tion. This survey has been built up by means of 580 
measurements of supernova magnitudes. To each magni¬ 
tude p. one associates the corresponding redshift z and 
magnitude error 5p, assuming the redshift error to be 
negligible. 

Under the hypothesis of a perfect Gaussian distribu¬ 
tion, it is possible to relate the luminosity distance Dl to 
the magnitude p. Afterwards, since Dl = Dl(Hq, qo,jo)-, 
one gets: 

Dl = —— (l + (IliZ + d,L2Z 2 + ■ ■ ■) , (14) 

no 

with the coefficients dm = i(l — qo) and dm = g(— 1 — 
jo + Qo + ^Qo) [3- Notice that we a priori assumed flk = 
0, without losing generality, since A does not depend upon 
it. 

We therefore consider a parameter-fitting procedure 
which makes use of a grid in which all parameters evolve, 
in their corresponding priors, i.e. Hq G [60,85], qo G 
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[—1,0], jo £ [—2,2], We minimize the chi-square func¬ 
tion: 

X 2 = A-^+V, (15) 

where A = x T C~ 1 x, B = JA(C _1 x)j, ^ = TrfC -1 ] 
and V = log(^r) and we get the results provided in 
Tab. [H The advantage of using the union 2.1 compilation 




FIG. 2: Contours for the experimental results. We plot the 
graphics of jo versus Ho (top left figure), qo versus Ho (top 
right), qo versus jo (button left) and jo versus Ho of the sec¬ 
ond fit (button right), in which qo is fixed in the framework of 
the ACDM model, with = 0.318, compatible with modern 
Planck results. In the contour plots, we indicate la (inside el¬ 
lipse) and 2a (outer ellipse) with different colors, correspond¬ 
ing to 68% and 95% confidence levels. 

on a grid is that previous systematics of other catalogs 
have been reduced. Thus, it is easy to show that the y- 
squared function (1151) permits to maximize the likelihood 
function: 

£ oc exp(-y 2 /2), (16) 

with a good approximation. 

In Tab. IU we present the results derived from our ex¬ 
perimental analysis. Fig. [2] contains the corresponding 
contour plots with the ler and 2<j of our coefficients. The 
values of qo and jo are basically close to each other, as 
theoretically predicted. Unfortunately, the higher errors 
on jo do not allow to characterize this property at a 5cr 
confidence level which would definitively suggest the va¬ 
lidity of our model. However, incoming refined data are 
expected to certify this experimental fact, to see how the 
coefficients qo and jo differ from each other. 


Parameter 

Fitl 

Fit2 

x 2 

0.975 

0.976 

Ho 

69.887+i : ™ 

69.873+?-*“ 

go 

—0.527+3°!! 

-0.523 === 

jo 

0-501+°;»5» 

0-482+O-558 

A 

—o.257+2 ; 8“ 

-0.255±g : g5| 


TABLE I: Table of our experimental results and numerical 
limits. We enumerate the coefficients qo and jo, showing that 
their absolute values are surprisingly close to each other, as 
predicted by our repulsion model in a FRW universe. We also 
evaluate A at present time, through the logarithmic propaga¬ 
tion of experimental errors. The subscript 0 indicates the cor¬ 
responding quantity is evaluated at present time. The fits are 
two: in the first one, we leave all coefficients free to vary, 
whereas in the second one we assume a fixed qo from the 
ACDM calibration, as suggested by Planck measurements. 


Conclusions and perspectives We presented a 
new approach to describe the universe dynamics at 
current and early phases of its evolution. We showed 
that it is possible to infer a process enabling the universe 
to accelerate by means of repulsive effects. To investi¬ 
gate such a phenomenon, we use an invariant definition 
of repulsive gravity in terms of the curvature tensor 
eigenvalues and we got two independent eigenvalues: A 
and t. It turns out that the dynamics of the repulsive 
component is governed by a particular eigenvalue A 
whose behavior at the moment of the repulsion onset, 
A = 0, allows us to construct a cosmological model 
based upon the behavior of the Hubble parameter at 
that moment. Moreover, the repulsion onset condition 
implies that the cosmographic parameters must satisfy 
the algebraic relationship j + q = 0, which is then 
shown to be in good agreement with observations. As 
a result of this analysis, we were able to feature an 
effective emergent dark energy fluid , which depends on 
two free integration constants whose values are fixed 
for dark energy and inflationary phases. In so doing, 
we demonstrated that this process fuels the universe 
acceleration without the need of scalar curvature flk ■ We 
therefore found a theoretical explanation for assuming a 
spatially flat universe at both late and early phases of its 
evolution. For these reasons, we believe our treatment 
to represent a third way to account the problem of 
cosmic speeding up, instead of considering barotropic 
dark energy fluids or modifications of general relativity. 
Finally, we checked the goodness of our cosmological 
approach by fitting the union 2.1 compilation, inferring 
cosmic bounds on the observable quantities and showing 
that our picture is compatible with present time obser¬ 
vations. Future developments will better characterize 
our model at different stages of the universe expansion 
history, to better clarify how eigenvalues change sign. 
Additional efforts will be also spent to better constrain 
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the cosmological parameters of interest with different 
data sets and combined data surveys. 
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